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We calculate ordinary and radiative muon capture on the proton in an effective field theory of 
pions, nucleons and delta isobars, working to third and second order in the small scale expansion 
respectively. Preceding calculations in chiral effective field theories only employed pion and nucleon 
' degrees of freedom and were not able to reproduce the photon spectrum in the pioneering experi- 

X^"^ , ment of radiative muon capture on the proton from TRIUMF. For the past few years it has been 

. speculated that the discrepancy between theory and experiment might be due to A(1232) related 

effects, which are only included via higher order contact interactions in the standard chiral approach. 
, In this report we demonstrate that this speculation does not hold true. We show that contrary to 

' expectations from naive dimensional analysis isobar effects on the photon spectrum and the total 

' rate in radiative muon capture are of the order of a few percent, consistent with earlier findings 

in a more phenomenological approach. We further demonstrate that both ordinary and radiative 
muon capture constitute systems with a very well behaved chiral expansion, both in standard chiral 
perturbation theory and in the small scale expansion, and present some new ideas that might be 
at the bottom of the still unresolved discrepancy between theory and experiment in radiative muon 
capture. Finally we comment upon the procedure employed by the TRIUMF group to extract new 
information from their radiative muon capture experiment on the pseudoscalar form factor of the 
nucleon. We show that it is inconsistent with the ordinary muon capture data. 
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I. INTRODUCTION 



With the ubiquitous success of the Standard Modelg ordinary and radiative muon capture on the proton can 
nowadays be considered as an excellent testing ground for our understanding of spontaneous and explicit chiral 
symmetry breaking in QCD. This stems from the fact that the typical momentum transfer in these reactions is very 
small — of the order of the muon mass — and one therefore can apply effective field theory methods, in particular baryon 
chiral perturbation theory. 

Ordinary muon capture (OMC), 

fi-{l)+p{r)^i^f,{n+n{r') , (1.1) 

where we have indicated the four-momenta of the various particles, allows to measure the so-called induced pseu- 
doscalar coupling constant, gp. This coupling constant is nothing but the value of the induced pseudoscalar form factor 
Gp{t) at the four-momentum transfer for muon capture by the proton at rest, gp = Gp{t = — 0.88m^)/2MAr, with 
n^fi (Mn) the muon (nucleon) mass. Theoretically, it is dominated by the pion pole as given by the time-honored 
PCAC prediction, g^'^^'^ — 8.89. Adler and Dothan |^ as well as Wolfenstein Q calculated the first correction to the 
PCAC result utilizing by now outdated (and in some cases misleading) methods. The ADW relation was rederived 
directly from QCD Ward identities within the framework of heavy baryon chiral perturbation theory |^ and shown 
not to be affected by A(1232) effects in [||, i?™^'^ = 8.44 ± 0.23 (a similar calculation with a slightly different result 
was given in 1^). The presently available data have, however, too large error bars to discriminate between the pion 
pole prediction and its corrected version ||^. Furthermore, Gr,{t) for low four-momentum transfer squared can be 
extracted from charged pion electroproduction measurements [^1 . The resulting momentum dependence of the induced 
pseudoscalar form factor is in agreement with the pion pole prediction. However, as in the case of the pseudoscalar 
coupling, the data are not precise enough to be sensitive to the small corrections found in |H [|| j|] . 
While the momentum transfer in OMC is fixed, radiative muon capture (RMC), 

fr{l)+p{r)^,.^{l') + n{r')+^{k) , (1.2) 

has a variable momentum transfer t and one can get up to t = at the maximum photon energy of about 
k ~ 100 MeV, which is quite close to the pion pole. This amounts approximately to a four times larger sensitivity 
to gp in RMC than OMC. However, this increased sensitivity is upset by the very small partial branching ratio in 
hydrogen (~ 10^^ for photons with fc > 60 MeV) and one thus has to deal with large backgrounds. Precisely for this 
reason only very recently a first measurement of RMC on the proton has been published ||9[ . The resulting number 
for gp, which was obtained using a relativistic tree model including the A-isobar | [Tc| ], came out significantly larger 
than expected from OMC, g^^'^ = 12.35 ± 0.88 ± 0.38 = 1.46 5^"^"^. It should be noted that in this model the 
momentum dependence in Gp{t) is solely given in terms of the pion pole and the induced pseudoscalar coupling is 
obtained as a multiplicative factor from direct comparison to the photon spectrum and the partial RMC branching 
ratio (for photon energies larger than 60 MeV) . It was also argued in ||^ that the atomic and molecular physics related 
to the binding of the muon in singlet and triplet atomic and ortho and para p/ip molecular states is sufficiently 
well under control. 

The TRIUMF result spurred a lot of theoretical activity. While radiative muon capture had already been calculated 
in phenomenological tree level models a long time ago, see e.g. (ll| [|l^ ll^ [|l^, heavy baryon chiral perturbation 
theory was also used at tree level including dimension two operators |]15| and to one loop order fl^ . The resulting 
photon spectra are not very different from the ones obtained in the phenomenological models, the most striking feature 
being the smallness of the chiral loops [ p^ , hinting towards a good convergence of the chiral expansion. At present, 
the puzzling result from the TRIUMF experiment remains unexplained. It is, however, a viable possibility that the 
discrepancy does not come from the strong interactions but rather is related to the distribution of the various spin 
states of the muonic atoms. It is therefore mandatory to sharpen the theoretical predictions for the strong as well as 
the non-strong physics entering the experimental analysis]^ 

Here, we wish to reanalyze RMC in the framework of the so-called small scale expansion pj|j20|], which allows 
to systematically include the A resonance into the effective field theory. Although Ando and Min [|16[ have already 



'^For a recent review of utilizing various muon capture reactions as filters for physics beyond the Standard Model see [Q]. 
recently proposed solution to the problem based on a novel term has been shown to be inconsistent with CYC in 
ref. ^|. The corrected term- which in the chiral counting comes in at C'(p'')-only gives a small contribution and is already 
contained in the RMC calculation of Ando and Min [hdl. 
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shown that the RMC process possesses a well behaved chiral expansion up to N^LO, it has been noted quite early |2|] 
that one should reanalyze RMC in a chiral effective field theory with explicit A degrees of freedom^. This is due to 
the fact that the A-resonance lies quite close to the nucleon and therefore, in a delta-free theory like HBChPT as used 
in , could lead to unnaturally large higher order contact interactions which would spoil the seemingly good chiral 
convergence. Stating the same concern in the language of (naive) dimensional analysis, it suggests the possibility of 
corrections of the order of 30% due to the small nucleon-delta mass splitting, mp/(AjfA — Mm) ^ Sm^/A/jv- In the 
small scale expansion, the leading delta effects involving the large Ml ^N/S. vertex already appear at second order 
(e denotes a genuine small parameter, the pion mass, external momentum or the A^A mass splitting) and can 
therefore already be analyzed at tree level. The study presented here therefore constitutes a natural extension of the 
two previous analyses of RMC using chiral effective field theorie^ Some preliminary results where reported 

in ref . p3[ | . As we will show, the small scale expansion allows for a very transparent separation of the resonance and 
chiral pion effects. 

The manuscript is organized as follows. In section 2, we briefiy review how the Standard Model at low energies is 
mapped onto a chiral effective field theory. The ingredients of this field theory, which uses pions, nucleons and the 
delta isobar as degrees of freedom, are discussed in section 3. Section 4 is concerned with ordinary and radiative muon 
capture on the proton. In the framework of the small scale expansion in OMC the leading order delta effects only 
come in at N^LO (i.e. 0{e^))^ whereas for RMC one can already study the explicit infiuence of this resonance at NLO 
(i.e. ©(e^)). The status of the induced pseudoscalar form factor is reviewed and the analysis of the TRIUMF RMC 
experiment is critically assessed in section 5. Section 6 contains the summary and conclusions. Some more technical 
aspects of this work are relegated to the appendix. 



II. FROM THE STANDARD MODEL TO THE EFFECTIVE FIELD THEORY 



In this section we briefly summarize how QCD coupled to the standard electroweak theory is mapped onto the 
pertinent effective field theory. For that, we start with the electroweak Lagrangian. For our purpose, we only need 
the coupling of the charged currents to the charged massive gauge bosons {W'^) and the coupling of the massless 
photon (A^) to the electromagnetic (em) current. 
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(2.1) 



where the last term refers to the non-leptonic weak interactions, gi and 52 are the gauge couplings of the U(l)y and 
the SU(2)i gauge groups and 0w is the weak mixing (Weinberg) angle. In terms of the light quark and lepton fields 
the em and charged weak currents read 
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Jch = ""'7^ (1 - 75) + v^,Y (1 - 75) M + 



(2.2) 
(2.3) 



with = {uw, dm) denoting the quark eigenstates of the weak interaction and the ellipsis denoting terms we do not 
need in what follows. For constructing the effective field theory, it is most convenient to consider the QCD Lagrangian 
coupled to these gauge fields treated as external local sources, which transform locally under chiral symmetry p4| . 
With e — gi cos 6w and = V^(a^) — -^aix), QCD coupled to these external sources takes the form 



q[s{x) ~ip{x)]q 



+q M°\x)- A^°Hx)l5 q-q s^°'> (x) - i p^°\x) q 



(2.4) 



^Phenomenological models have claimed for a long time that the A contribution does not exceed 8% in the photon spectrum 
for photon energies above 60 MeV [Q. 

*At the moment there exists no systematic scheme in chiral effective field theories to simultaneously include the effects of 
explicit vector mesons and other strong short range effects into the calculation. Given the small momentum transfer of RMC 
in the t-channel one expects, however, that these contributions can be correctly accounted for via Ofp^) contact interactions 
[jf6|. A Born term analysis of vector meson contributions in the RMC process has been presented in p2[. 
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with q the bi-spinor of the hght quark fields in the strong interaction basis. For the case of RMC we are confronted 
with the following scenario of external source^: 



s{x) = 
sW(a;) = 



, 

rhi , 



p(a;)-p(°)(a;) = 0, 



giVud ( 1 
V8 



V. {x) + h.c. 



.g2Kd / 1 1 



V8 



(x) + h.c. 



= , 



(2.5) 
(2.6) 
(2.7) 

(2.8) 
(2.9) 

(2.10) 
(2.11) 



with s(x)^p{x)^Vy,(x) and ap,(x) scalar, pseudoscalar, vector and axial-vector fields, in order, in an obvious isospin 
(singlet and triplet) notation. V^d is the pertinent element of the CKM matrix. Note that the term >Cqcd Eg. (p^) 
is chirally symmetric. The explicit chiral symmetry breaking due to the current quark masses resides in the zeroth 
component of the external scalar source. 

Having specified the external field environment we now analyze the required strong matrix elements for calculating 
OMC and RMC. First we define quark vector and axial- vector currents 



(2.12) 
(2.13) 



Here, the T"^ — t°- /2 are the generators of SU(2). With these definitions one then specifies the strong matrix elements 
which have to be calculated in the effective theory 



a) {N\V;\N) , 

b) {N\Al\N) , 

c) {N\TV;^V^\N) , 

d) {N\TV"^Al\N) , 



(2.14) 
(2.15) 
(2.16) 
(2.17) 



with all possible combinations of the external fields of Eqs.(p^- 2.11 ) and T denotes the conventional time-ordering 
operator. To proceed, we now have to specify the effective Lagrangian which will be used to calculate these matrix 
elements. 



III. CHIRAL LAGRANGIANS 

In this section, we briefly review the chiral Lagrangian underlying our calculation. While the meson part is standard, 
with the external momenta and the pion (light quark) mass counted as small parameters, for the meson-baryon system 
we include the nucleons as well as the A(1232) resonance. To systematically account for the effects of the latter, the 
iVA mass splitting is considered as an additional small parameter. This is routed in phenomenology and the large Nc 
limit of QCD, but not in the strict chiral limit, where the A decouples. These three small parameters are collectively 
denoted by e. 

A. Meson chiral perturbation theory 

At low energies the chiral symmetry of QCD is spontaneously broken to the vectorial subgroup SU{2)v'- SU{2)l x 
SU{2)i^ — !■ SU(2)v- The strictures of the spontaneous and the explicit chiral symmetry breaking can be explored in 



'We are working in the limit of no isospin breaking, i.e. equal quark masses niu = rud and no internal (virtual) photon effects. 
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terms of an effective field theory, chiral perturbation theory (ChPT). As a tool, one works with an effective Lagrangian, 
which consists of a string of terms with increasing dimension. For our purpose, we only need the first term in this 
expansion, the non-linear cr-model chirally coupled to the external sources. 



(2) 

TTTT 



with 



U{x) 



exp <; — T 



t:{x) 



X = 2Bo s + s 



= (0) 



ip + ip 



(0)^ 



(3.1) 

(3.2) 

(3.3) 
(3.4) 



The pions, which are nothing but integration variables, are collected in the matrix-valued field U{x) and the explicit 
symmetry breaking due to the light quark masses is hidden in x via the zeroth component of the scalar source. Fq is 
the (weak) pion decay constant (in the chiral limit) and Bq is related to the scalar quark condensate. We work in the 
standard scenario with Bq ^ F^. This specifies completely the meson part of the effective Lagrangian. 

B. Including baryons: The small scale expansion 

The nucleon-delta-pion system chirally coupled to the external fields can also be represented by a Lagrangian, 
which decomposes into a string of terms with increasing dimension. We work here in the heavy mass formulation, in 
which the nucleon and the delta are essentially considered as heavy, static sources. This allows to shuffle the baryon 
mass (Mb) dependence into a string of 1/Mb suppressed vertices and gives rise to a consistent power counting p^ ] 
p6[ [ pO| . Denoting by N the (heavy) nucleon isodoublet field and by the Rarita-Schwinger representation of the 
heavy spin-3/2 field, the lowest order terms read 

'(1) 



C);>,^N[tvD + gAS-u]N 



(1) _ 

ttA 



ft[tv D'^ - S'^Ao + ...] g^, , 



^^^NA = 9-NA {ft g^a wfN + N W. 



T! 



(3.5) 



with Aq = Ma — Mq (Mq) being the nucleon-delta mass splitting (nucleon mass) in the chiral limit (to the order we 
are working, we can set Aq = A). All other quantities Q in the chiral limit are denoted as Q. Furthermore, and 
5^ are the four-velocity and the spin-vector of the heavy nucleon. The various chiral covariant derivatives, chiral 
connections and vielbeins are (for details, we refer to and pfl]) 















Mm 













(3.6) 



where i, j, k are isospin indices. At second order in the small scale expansion, we need the following terms 



(2) 

'ttN 



2M, 



,(2) _ 
-ttNA — 



■N{{vD) 



igA (S-Dv-u + v- uS-D) 



2Mn 



[bi ^f+^,S" + . . .] TV + h.c. 



N 



(3.7) 



(3.8) 
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with 



F^u = d^F^ ~ d,F^ ~ I [F^^F^] ■ X = L,R, 
FH = V,. + , Fl: = - a^ , 

v^^} = dA'^ ~ d.v\:^ . (3.9) 

Note that in C)^j^ only two low-energy constants (LECs) appear, which we have expressed in terms of the isoscalar 
and isovector anomalous magnetic moments, Kg and k^,, of the nucleon, respectively. The LEG hi is fixed from neutral 
pion photoproduction at threshold, hi = 12.0. The C(e^) Lagrangians given so far are sufficient to calculate the 
leading order A(1232) effects in RMC. However, as mentioned above, in OMC the leading A(1232) related effects 
only occur at N^LO — specifically as a C(e^) contribution to the form factors of the nucleon. In SSE these terms have 
already been calculated in Q and we refer the reader interested in the details of the third order Lagrangian required 
for such a calculation to ref . Q . We now have specified the effective Lagrangian needed to work out OMC and RMC 
to second order in the small scale expansion. The pertinent Feynman rules to perform these calculations are collected 
in appendix A. 

IV. MUON CAPTURE 

This section is concerned with the theoretical description of OMC and RMC. To keep the manuscript self-contained, 
we give all formulae necessary to calculate the capture rates. 

A. Leptonic matrix elements 

Consider first the purely leptonic part of the muon capture reaction. For OMC and RMC, we need the following 
two leptonic matrix elements: 
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{^^\J+\^i) = ^^^p,{nl,i^-l5)K^) , (4.1) 



(^M 7l J+Ia^> ^-'^^ Ml'h^. (1 - 75) (2e* • /- ji*) , (4.2) 

with e* t he p olarization vector of the photon. The corresponding four-momenta have already been defined in Eq.( |l.l| ) 
and Eq.(1.2), respectively. 



B. Hadronic matrix elements 



For OMC one only needs to calculate the vector and axial-vector transition matrix elements (7i|V \p), (nlA \p) of 



Eqs.(2.14 2.15). In the limit of exact isospin symmetry these matrix-elements constitute the off-diagonal components 
of the isovector nucleon vector and axial- vector currents which are typically parameterized in terms of four form factors 
(e.g. ref. ^): the isovector vector form factors F^2i '^^11 as the axial and the induced pseudoscalar form factor. 



GA{q^) and Gp{q^) respectively. In terms of these, the relativistic strong matrix elements defined in Eqs.( 2.14 , 2.15 ) 
(which corresponds to the vector and axial correlators, shown in the first row of fig.p]) are given by: 



{n\V^ \p) = -i^^^n{p2) 



F^iq'h, + ^FUq')<y,.q'' 



GA{q^)l^^b 



Gp{q' 



2M 



'J-N 



-q^iib 



p{pi) 
p{pi)- 



(4.3) 
(4.4) 



where q^ = (p2 — Pi)'^ is the invariant momentum transfer squared. The Dirac F^^^ and Pauli i^2"'' form factor are 
subject to the normalizations 



FtH0) = l, Fr{0)^n., 



(4.5) 
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with K„ = 3.71 the isovector nucleon anomalous magnetic moment. In the axial matrix element, Eq.(4.4) one has 
assumed the absence of second class currents. The electromagnetic form factor are rather well known. The current 
situation concerning their theoretical understanding and experimental knowledge can be found e.g. in Muon 
captures therefore provide us with the opportunity to study the weak axial structure of a nucleon. While GAiq^) can 
be extracted from (anti)neutrino-proton scattering or charged pion electroproduction data, Gp{q^) is harder to pin 
down and in fact constitutes the least known nucleon form factor. In Fig.|^ we present the "world data" for Gp{q^). 
These four form factors have already been calculated to 0{e^) in One findsQ 



i+g(r-r)V + o(g') 



l + -{rA)V+0{q') 
6 



2n -iMNgTTNNFTr ^ 71 #2 2 
'^Pil ) = — — 12 gf^Mjvr^ 



^2 - gf2 



(4.6) 
(4.7) 
(4.8) 



where we have systematically shifted all appearing quantities to their physical values. 



' 1.2 



and rA are the isovector 



Dirac and Pauli radius and the axial radius respectively. They receive contributions from chiral loops and counter- 
terms, except for which is free of any LEG. Detailed results can be found in In Fig.^ the difference between 



the usual pion-pole parameterization for Gp(q^) and the full chiral structure of the form factor, Eq. ^.8| is displayed. 
Note that this chiral structure is not affected by A(1232). In the kinematical region of RMC, which mainly lies to 
the "left" of the OMC point in FigJ^ the structure effect in Gp{q^) proportional to the axial radius ta is expected 
to play only a small role. Certainly, the present experimental uncertainties both in OMC and in RMC 1^ are too 
large to distinguish between the two curves, but new ef forts are under way 

Usually the non relativistic reduction of Eqs.(4^,4^) is done in the Breit frame. Here we give the non relativistic 
strong matrix elements in the rest frame of the proton where our calculation is done: 



(^iVTb) 



{n\A-\p) 



g2Vud .r - / N 

-i — — A/2 ny{p2) 



^M^Fi^\q')-^^F^^\q^)\v, 



E2 + IVl N 

2 



V8 i\E2 + Mn' ' ' E2 + MN 

Lj-(^,->(,^)+f-->(,^))-^fr(,',],„ 

|S„,S-5|(f<"'(,») + F«(,»))}p„(0), 

2 S* • (jwp 



E2 + MN 

. g2Vud 



Af2ny{p2)<GAiqn 



+RGpiq^) 



S ■ qq^ 



Mn {E2+Mn) 



2 5. - 

^ E2 + Mn 



(4.9) 



(4.10) 



where J\f2 is the usual normalization factor of the neutron wave function, A/2 = y ^2Ai^" ^^^^ ^2 is the neutron 

energy. One gets immediately the results to 0{e^) by replacing the form factors by their values at = 

We now turn to the vector-vector (VV) and vector-axial (VA) correlator which we only need to C(e^) in SSE. 
Working in the Coulomb gauge e* ■ v — for the photon and making use of the transversality condition e* ■ k = 0, we 
find (the pertinent Feynman diagrams for the VV correlator are shown in the second row in fig. |l] and the ones for 
VA in figs. II) 



^Note that due to the small momentum transfer < 0.01 GeV^ in OMC/RMC it is sufficient to work in the radius 
approximation of the form factors, i.e. to truncate the g^-dependence after the first term. The full 0(€"*) momentum dependence 
for -q^ < 0.2 GeV^ can be found in 1). 

^ Let us briefly emphasize that there exists another window on Gp{q^) — pion electroproduction. So far there has only 
been one experiment that took up the challenge, with the results shown in Figj^. In this kinematical regime the structure 
proportional to rA produces the biggest effect and a new dedicated experiment should be able to identify it — thereby enhancing 
our knowledge of this poorly known form factor considerably! In fact, at the Mainz Microtron MAMI-B a dedicated experiment 
has been proposed to measure the axial and the induced pseudoscalar form factors by means of charged pion electroproduction 
at low momentum transfer Ipsll. 
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\TV-e*V-\p) 



{n\TV -t^A-Jp) 



(2) 



^— n^yr ) 



1 



1 



IN 



MnU! 



(2) ■ 92yud e _ / /X , 

^ ' = —I ;= — n„(r ) X 



2RgA 



S-{r' - r) 

(r' — r)2 — 



[(1 + Ki,) [S* • £*, • k] 



2e*-il-l') (l-V), 



r]+0(l/M^) p„(r) 



P 9A {v ■ r' - V ■ r) S ■ {r + r') 



{I - I'Y - ml ^_ 
2e* • 



M 



N 



-2RgA 

9A 

~ Mn 



(r' - r)2 
vl 



mi: 



2Mn 



{I -I') 



t\2 



mi 



9a „ * 



(2 + + At^) 5'" [S-t*,S-k] ^ {K^-K,)[S-e\S-k] 



2 5"e* -r 



+ 



9itN Abi 



3M. 



N 



2A[k°'S-e 



9fia — R 



(l - I'Y - ml 



ujv"S -e* ~ e*°'S ■ k] 4 S"[S ■ e* , S ■ k] 



A2 - 

4[S-e*,S-k] g° 
3(A- w) 



3(A + w) 



(4.11) 



(4.12) 



with uj = v-k and i? = 1 in QCD. We have introduced this factor multiplying the Born term contributions proportional 
to the induced pseudoscalar form factor for the later discussion. One can easily check from the continuity equations 
satisfied by the correlators (which for example relates the vector-axial correlator to the axial one, sec |22|) that gauge 
invariance is satisfied in the above equations. We note again that here we only give the results to O(e^) as this 
is sufficient to study the leading order A(1232) effects in RMC. We also want to point out that the vector-vector 
correlator is free of delta effects to this order 0{e^), i.e. the leading A(1232) effect only appears in the vector-axial 
correlator, cf. fig ||. 



C. Ordinary Muon Capture 



Surprisingly, no full analysis of OMC exists in the literature for the framework of chiral effective field theories. 
All previous analyses 1^-^ stopped at the level of writing down the relevant matrix-elements/form factors (here 
Eqs.(4.3 4^)), but no discussion of the implications for the lifetime of OMC was given, which can be written down in 
a closed analytic form. In this section we are going to fill this gap. We work in the Fermi approximation of a static 
field, i.e. the gauge boson propagator is reduced to a point interaction (since the typical momenta involved are 
much smaller than the W mass), 



M 



OMC 



* TIT [("IKTb) - {n\A-\p)] 



(4.13) 



1. Spin-averaged OMC 



Introducing the Fermi constant Gf via Gf = .92 V^/(8Af^), we define the square of the spin-averaged invariant 
matrix element to be 



OMC|2 _ '■^F^ud T(a) TrfU' 



(4.14) 



With the normalizations 



7 



■■^ — ' 2 m^- ^ — ' 



5]p„(r,a)p.(r,a)=i'+(l + |^j, ^ a)n„(r', a) = F+ h + |^ ) , (4.15) 



one then obtains the symmetric tensors 



L^^} = {1^.11 - gt.J ■ I' + I'^l. + ie^c..pl"l"^} , (4.16) 



Hl^J = a^i^^u + ^ {-"p-^v - + ct'/^ + d (^;^r:, + r'^v,) + eiefy^vp . (4.17) 

Note that we have evaluated the tensors for the special kinematic condition of both the proton and the muon being 
at rest^, i.e. = (m^-, 0,0,0), — (0,0,0,0). a,b,c,d,e are coefficients which depend on the four form factors 



Eqs.(4.6, 4.7). For illustration we give them to 0{e ) 



We now assume that the initial muon-proton system constitutes the ground-state of a bound system. We therefore 
replace the plane-wave wavefunction of the muon used so far in the calculation by the Is Bohr-wavefunction ^(x)is of 
a muonic atom. We note that to 0{e^) we are not sensitive to the extended structure of the proton - any dependence 
of the capture process on the electric/magnetic radius of the nucleon will only occur at 0{e'^) As we will see 
below, the terms of O(e^) in OMC only represent a small N^LO correction. For simplicity we therefore approximate 
the muonic atom wave- function dependence by its value at the origin even at 0(6^), which effectively constitutes an 
upper bound for the 0{e^) correction to the total width. In a systematic 0{€^) analysis one would of course have 
to calculate the overlap between the Bohr-wavefunction and the respective proton-neutron transition form factors. 
Confining ourselves to a; = we use 

Q,3/2„3/2 

$(0)1, = ^ , (4.19) 

with the reduced mass /i — AlNm^- / (Mn + ^7ip-) and a = e^/ATrhc in the Heaviside-Lorentz convention. We can 
therefore calculate the spin-averaged rate of ordinary muon capture via 

aa' ss' 

with the normalization factors Ju — 2Ei,, Jn = 1 + v -r' /Mn ■ Evaluating all expressions at 0{e'^) accuracy one obtains 

ToMC = (247.0 - 61.6 - 4.0 + 0{1/M}^)) x s'^ 

= 181.5 X s-^ (4.21) 

Here we have used the coupling constants, masses and radii as given in table ^. r^"' and va corresponds to their 
empirical values since the value of the counter-terms which enter in their expressions to O(e^) has been set to 

reproduce them. The value of (?'2'''')sse ~ ^-^^ given counter-term free by the 0(6^) SSE calculation. For 

comparison we also give the spin averaged OMC capture rate to 0{p^) in HBCHPT without explicit spin 3/2 degrees 
of freedom: 



ToMc'''^^ = (247.0 - 61.6 - 3.8 ©(l/M^)) x s^^ 

= 181.7 xs~^ (4.22) 



This approximation is well-justified due to the low binding energy E « 2.5 keV of an s-wave muonic atom as compared to 
the muon mass. 
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The only change occurs at N^LO via (?'2'''')HBChPT ~ ^-^^ fm^, showing only a very weak dependence on the exact 
value of the isovector Pauli radius, which has the "exact" value ir^^)l^p — 0.80 fm^. We note that in Eqs.(4.21) 



the 0{e^) SSE contribution amounts to a correction of 25% of the leading term, while the 0{e^) one is two orders 
of magnitude smaller, indicating that muon capture on a proton really constitutes a system with an extremely well 
behaved chiral expansion. 

Finally we discuss the case of no explicit chiral symmetry breaking (i.e. in the chiral limit m,r = 0). One expects 
the spin-averaged capture rate to behave as 

r^MC = (214 - 46 +0(1/A4)) X s-i 

= 168 X s"i . (4.23) 

which only represents a small finite shift compared to the value for finite quark masses. The leading infrared singu- 
larities only occur at N^LO due to the long range nature of the chiral pion cloud in the isovector form factor radii 

i- 

The reason for the nice stability of perturbative calculations for OMC in the physical world of small finite quark 
masses is of course the fact that contributions of order n are suppressed by {rrii/K^)^ ^, with i — Tr,fj, and A,^ ~ 
Afjv ^ 1 GeV. Note that one could also retain the higher order kinematical corrections starting at order 1/M^. To 
be more precise, this refers to the energy-momentum relation between the various particles, the mass term appearing 
in the various projection operators and the phase space. In that case, only the various correlators (A,V,AV,VV) are 
truncated at ord er l/A /jy. Using this approximation, the total rate is Fomc = 193 s^^, not much different from the 



one given in Eq.(4.21). We remark already at this point that in the case of RMC, it is mandatory to retain these 



higher order terms if one stays at a low order in the small scale expansion as done here. 



2. Hyperfine effect in OMC 

The spin-averaged OMC scenario presented in the previous section is only of theoretical interest. In nature the weak- 
interactions shows a very strong spin-dependence, which leads to quite different decay-rates depending on whether 
the captured Is muon forms a singlet or a triplet spin-state with its proton pT[ |, where the singlet is the usual state 
(l/-\/2)(| tii) — I lA)) in terms of the muon and proton spins and the triplet accordingly. Although the hyperfine 
splitting between the two levels "only" amounts to 0.04 eV, the occupation numbers of the levels due to thermal and 
collision induced processes tend to be far from statistical equilibrium. In order to make any contact with experiment, 
the singlet/triplet rates need to be calculated separately. These various spin states are obtained by using the following 
projection operators for the muon, 

MU/2MU/2) = i(l+75.?)''^+™" 



2' ' 2m^ 



/i(Z,±l/2)/2(/,Tl/2) = — --,5707:ili±il2hoil/+m^) , (4.24) 



and similarly for the proton, 

uil, l/2)u{l, 1/2) = 1(1 + 75 ^ 1(1+ V) (^1 + , 



u{l,±l/2)u{l,Tl/2) = 1(1+ V) 2Mn{eT+Mn) ^^^4^^^ ^ (y^+2M^) 1(1+ ^ . (4.25) 

For the total capture rates of singlet and the triplet states in the muonic atom, we then find the following decomposition 
into leading, next-to-leading order pieces and next -to-next-to-leading order pieces to 0{e^) in SSEjj: 

F^'j^fc = (957- 245GeV/A/Ar + (30.4GeV/M^ - 43.17) + 0{l/M^)) x s^^ = 687.4 x s^^ , 

FojPp = (10.3 + 4.72 GeV/MAT - (1.22 GeV/M^ + 1.00) + ©(l/M^)) x s"^ = 12.9 x s"^ , (4.26) 



®For the remaining part of the discussion on OMC we only give the values calculated in to 0{e^) in the small scale expansion. 
The corresponding values to 0{p'^) in HBChPT are very similar. For OMC the two effective field theories only differ in the 
N^LO contribution of the isovector Pauli radius as discussed in the previous section. 
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displaying the dramatic spin- dependence due to the V-A structure of the weak interaction in the Standard M ode l 
These numbers correspond to a value of g-nNN as given in table | and includes the pion pole corrections, Eq.(4 



Since Gp contributes negatively to the singlet rate a larger value of g-^NN leads to a smaller value for the rate: 
rp MH = 681.9 X s~^ for g^^NN = 13.4. Similarly, neglecting the pion pole corrections leads to Fq^^^ = 676.1 x s~^. In 



Eq.(4.26) we have split the 0{e^) term (third and fourth terms in parenthesis) into the contribution from the 1/M^ 
corrections and the pure C(e'^) terms stemming from the various radii which lead to the dependence of the form 



factors, Eq.(4.8). It is very interesting to note that these two contributions more or less cancels themselves in the case 
of the singlet term. It is thus extremely important to perform a consistent chiral expansion. In a relativistic Born 
model one obtains the simple nice formula^ for Tq^^ |31] : 

ToMC - (6.236 Fi"(g2) + 0.5513 ^^^"(g^) + 16.44 G^(g2) _ 0.2834 Gp(g2))' (4.27) 
- 683 X s~i (4.28) 

where q^ = — O.SSm^ is the momentum transfer and the remaining parameters are taken from table ^. The smaller 
value quoted in |30|, ||3l| comes from the fact that in the sixties gA was somewhat smaller and rA somewhat bigger. 
This formula though very appealing should not be used, being in contradiction with the modern viewpoint of power 
counting. The good agreement with the SSE result is purely accidental. The problem arises from the fact that Tq^^^ 



is a rather sensitive quantity as can be seen for example in Eq.(4.27). Indeed the terms proportional to -^2(90) ''^^'^ 
Gpilo) ^^'^ of same order of magnitude but have different signs so they have a tendency to cancel each other 
rendering the values of F'^^^s rather sensitive to the exact values of these two quantities. 

So far we have only considered OMC for the case of muonic atoms. For the case of a liquid hydrogen target one 
also has to take into account the possibility of muon capture in a muonic molecule pfip, which can be formed via the 
reaction + pep — ^ pfip + e + 124eV. In such a molecule the muon can be found in a so called ortho (O) (spin of the 
protons parallel) or para (P) (spin of the protons antiparallel) spin state relative to its two accompanying protons. 
The decay rates of these molecular states can be calculated from the singlet /triplet rates of the muonic atom via 



1 

■4' 

To = J3l/2ri/2 +P3/2r3/2 , (4.30) 



Tp = 27P-(3rt,ip + r,i„g) , (4.29) 



The wavefunction corrections]^ are taken to be 70 = 0.500, 7p — 0.5733 We note that the para molecular state 
is often referred to as the statistical mixture, as it corresponds to the naively expected occupation numbers of the 
muonic atom. For a precise calculation of the ortho molecular state on the other hand one first has to know the exact 
probabilities Pi/2,P3/2 for the muonic molecule being in a total spin S=(l-F/2)=l/2 or a total spin S=(l+l/2)=3/2 
state, with > 0.5 |32j. The corresponding decay rates are given by ]3C| , |3^ 



ri/2 = 270 [^^Fsing + ^Ftrip j 

F3/2 = 270 FtHp (4.31) 



Theoretical calculations of the spin-effects in the muonic molecule |33 suggest P1/2 ~ 1, P3/2 ~ which leads to the 
values: 

fO^c ^ 208 X s-i 

T'^^'^ = 493 ■ ■ ■ 519 X s-^ . (4.32) 

where the range given in Fq'^*-" corresponds to 0.95 < Pi/2 < 1. We have allowed here for a 5% uncertainty in the 
occupation numberg|^ to show the sensitivity of our results on this quantity. Since Fging ^ Ftrip for OMC, Fq^^ 



^"We note that we have different relative signs compared to ref. (glj to translate into our notation. 

^^^P (70) denotes the ratio of the probability of finding the negative muon at the point occupied by a proton in the para- 
muonic (ortho-muonic) molecule and the probability of finding the negative muon at the origin in the muonic atom. We are 
grateful to Shung-Ichi Ando for pointing out ref. |28[ to us. 

^^The possibility of very different occupation numbers has been raised by Shung-Ichi Ando during the Chiral Dynamics 2000 
conference. This will be reported in a forthcoming paper, see ref. |37[. 
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turns out to be roughly proportional to Pi/2- We note that our number for capture from the molecular ortho state 
agrees very well with the most recent measurement T'^^' = (531 ± 33) x s~^ from Bardin et al. Let us stress at 
this point the importance of the new proposed experiment at PSI |Q which will be done with a hydrogen gas target 
and will thus be independant on these occupation numbers. One will directly measure Tsmg- 

With these expressions given above one can now calculate the rate for muon capture in liquid hydrogen for a general 
scenario, if one knows the relative occupation numbers for the atomic singlelj^ fs, the molecular ortho fo and the 
molecular para fp state: 

I'oMC ^ fs Tsing + /o To + fp Tp (4.33) 

For example, in the case of the recent TRIUMF experiment with fs = 0.061, fo — 0.854, fp — 0.085 one would 
obtain a total capture rate of FJ^^^^f ^ (504 + 0{1/Mfj)) x s^^. 

D. Radiative muon capture 

1. Total capture rates 

In the static approximation for the W-boson, the pertinent matrix element for RMC decomposes into two terms, 

A^M-p-^.n^ = {'^,\W+\^i)^^[{n\TV■e*V-\p) - {n\T V ■ e* A' \p)] 

+ {'^^l\W+\fi)i^[{nK'\p) - {n\A-\p)] , (4.34) 

so that its square in the spin-averaged case can be written as a sum of four terms, with both photons coming cither 
from the hadronic or the leptonic side and two mixed terms, i.e. 



aa' ss' W 



2r,2 y2 
RMC 1 2 _ ^ '^F^ud 



J (a) JTlll' , / T (b) TTt^iy , t {c) TTllU \ j {d) TTt^l' 



, AA' 



(4.35) 



with A, A' the photon helicities. Explicit expressions for the various tensors are not given here because they are lengthy 
and not illuminating. We also note that standard packages like REDUCE can not be used straightforwardly to obtain 
these tensors since the cyclicity of the trace in the presence of 75 matrices is not fulfilled. The total decay rate is 
given by: 

■'^("^^-'^ f\inedO pdc../-fl-f---(^--^mi \M---\\ (4.36) 



16^4 Jo Jo V V Mn JJ^ 



(7(7' SS 



AA' 



with Lo = ko the photon energy. The direction of the photon defines the z-direction and 9 in Eq.( 4.3(^ ) is the polar 



angle of the outgoing lepton with respect to this direction. The maximal photon energy is given by 



Furthermore, the energy of the outgoing lepton follows from energy conservation, 

= - ^ - ^ + ^ -^J^ i + 0{llMl) . (4.38) 

First we discuss the (academic) spin-averaged RMC scenario, which allows for a comparison with previous calculations: 



^^We assume that all muonic atoms initially in the triplet state are effectively converted into the atomic singlet state through 
collision with hydrogen molecules in the stopping target via the Gershtein-Zeldovich mechanism |^^. 
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rfpfnav. = (66.0+ 18.7 + ©(l/M^,)) X 10-3 =84.7x 10~^ (HBChPT) 

rfpfnav. = (66.0 + 20.4 + 0(1/^/1,)) X 10^3 =86.4x lO^^s"^ (SSE). (4.39) 

Both the HBChPT and the SSE results suggest a good convergence for the chiral expansion, as expected from 
dimensional analysis. Note that the leading order capture rates in both calculations are identical, as A(1232) related 
effects only start at sub-leading order. Our leading order result also agrees]^ with the HBChPT calculation of ref. 

However, our HBChPT 0{p^) correction is nearly 30% larger than the one given in [|15|. Comparing with the 
corresponding C(e^) correction in SSE, we note that A(1232) indeed leads to a larger decay RMC decay rate and 
constitutes a 9% correction to our 0{p^) contribution. However, the total (spin-averaged) decay rate is only affected 
by 2% due to the fast convergence of the chiral series for RMC. 

For the case of muonic atoms we obtain the following decay rates in the singlet/triplet channel|^, utilizing the 



projection formalism outlined in sec. [V C 



r^^^g^ = (12.7- 18.7 GeV/MAT + 0(1/71/^)) x lO'^s-^ = 3.10 x lO'^s'^ (HBChPT) 

r^^'^ = (12.7- 18.3 GcV/Mat + C'(1/M^)) x lO^^s^^ = 2.90 x lO^^s^^ (SSE) (4.40) 

r^RMC ^ (119- 3.86GeV/Afjv + C'(l/Af^)) x 10"^ s"^ 112 x lO^^s^^ (HBChPT) 

r^;^c ^ (119- 1.80 GeV/AfAr + C'(l/M^)) x IQ-^ s'^ = 114 x lO^^s^i (SSE) (4.41) 

Note that for the total numbers given we did not expand all kinematical factors in powers of 1/Mn since in case 
of the small singlet, the contribution from the terms starting at order 1/M^ can not be neglected. In fact, a strict 
truncation at 1/Afjv leads to an unphysical negative singlet capture rate. For the much bigger triplet, these higher 
order corrections are much less important. We remark that in ref. no strict l/M^ expansion was performed, only 
at some places these authors used the leading order results, e.g. for the nucleon energy by neglecting the recoil term. 
Only if one goes to a sufficiently high order in the small scale expansion, the truncation of these kinematical factors 
can be justified. Comparing the 0{p^) HBChPT with the C(e^) SSE calculation, we observe that while the total 
singlet capture rate is nearly identical in both approaches, the absolute value of the 1/Mm term in the total triplet 
capture rate is a factor of two different between SSE and HBChPT. However this NLO term is much smaller than the 
LO triplet rate, leading to a rather similar total spin triplet RMC rate with or without explicit A(1232) contributions. 

Finally we address the complications for RMC due to the presence of muonic molecules in the liquid hydrogen 
target. According to Eq.( 4.29| ), we can easily determine the capture rate from the molecular para state 



pRMC ^ 2 X IQ-^ s-^ (HBChPT) 

pRMC ^ gg_4 X g-i ^gsg<) ^4 42) 

Let us now turn to the molecular ortho state, which turns out to dominate in the recent RMC experiment from 
TRIUMF §. One obtains for pi/2 = 1: 

pRMC ^ 3Q 4 ^ g-i (HBChPT) 

Due to the triplet dominance in RMC (as opposed to the sing let dominance in OMC) Pg^^^ is now roughly proportional 
to (1 — 3/4pi/2) which leads to a big sensitivity of the RMC capture rate to the exact occupation numbers of the 
relative molecular sub-states. For example, a 5% uncertainty in the occupation numbers pi/2 = 0.95, — 0.05 
would lead to a 13% change in the ortho capture rate Fq^*-^ '--^ 35 x 10"^ . We will discuss the implications of this 
uncertainty when we compare our results with the measured photon spectrum from TRIUMF in the next section. For 
the total capture rate in the TRIUMF experiment 

^ RMC - /S 1 sing + /O I O + /P I P l4-44) 



^*Ref. fill gives rRMCHBChPT ^ ^ 10"^ s"^ to leading order. The small difference can be traced back to the different 

L^^^jj ^ spin 3j V . ^ 

values of some of their input parameters. 
^^Note the reversal of the relative size of the singlet to triplet contribution as compared to the case of OMC. 
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with fs = 0.061, fo = 0.854, fp = 0.085 § one would obtain T™(V'^'^ = (34.3 [34.8] + Oil/M^)) x lO'^s"! in 
HBChPT [SSE]. This leads to a relative branching ratio = FRMc/roMC 

^HBChPT ^ 34.3 X 10-3 ^ ^ ^ 0(l/i\4)) 
' 504 s ^ 

= ^'•\o4!"i'"' - X 10"' + '^(VM^)) . (4.45) 

Unfortunately the full relative branching ratio is not accessible in experiment, as one has to use a severe cut on the 
photon energies due to strong backgrounds. In the TRIUMF experiment only photons with an energy uj > 60MeV 
were detected. We therefore now move on to a discussion on the photon spectrum. 



2. Photon spectrum 



The photon spectrum dT/duj can be obtained straightforwardly from Eq.(4.36). We refrain from giving the lengthy 
formulae for the various atomic states here. We have calculated the photon spectra to 0{p^) in HBChPT and to 
0{e^) in SSE. The resulting curves are very similar. In figj^ we show the SSE 0{e^) results with the coupling values 
g-wNA X 6i = 1.05 X 12 = 12.6 for the singlet, triplet, para and ortho states. The relative difference between the 0{p^) 
HBChPT and the C(e^) SSE calculation for all states are shown in fig||, showing explicitly the small role of A(1232) 
in RMC. With the exception of the small singlet, the A effects amount to less than 5% for all photon energies. Only in 
the case of the singlet, a more pronounced influence of the A is observed. We note that these results are very similar to 
the ones found by Beder and Fearing (ijjl^ for the spectra and the relative contribution from the spin-3/2 resonance, 
although their calculation is based on a very different approach. Even the result for the singlet is comparable though 
not identical to the one of Beder and Fearing. It changes sign for photon energies of about 74 MeV. Like for the case 
of the small triplet in OMC, it is expected that the small singlet (for RMC) is more sensitive to 1/Mm corrections. 
We have also increased the coupling hi to values of 24 and 60 thus enhancing the A(1232) contributions in the SSE 
calculation by a factor of 2 and 5, respectively. We find very little sensitivity to this, e.g. the maximum in the photon 
spectrum for the para state changes from 1.47 GeV"^ s"^ for hi = 12 to 1.55 GeV"^ for bi = 60. Correspondingly, 
the total decay rate for the para state is increased by 7%. The smallness of the A(1232) related effects has two reasons. 
First, as noted already, it only appears at NLO, whereas the RMC process is mainly controled by the leading order 
effects, hinting to a good convergence behavior of the chiral expansions. Second, its sole contribution at that order 
is in axial-vector correlator, but not in the other three correlators. This already leads to a statistical suppression 
as compared to the pure nucleonic contributions. Third, despite the largeness of the coupling gT^NA x hi (remember 



that bi is related to the dominant Ml NA transition |36|), the delta contribution to the VA-correlator is still smaller 
than the one from the nucleon, which is enhanced by the large isovector magnetic moment. The dominant axial and 
axial-vector contribution comes indeed from the leading pion pole, which is not affected by the A(1232) effects. 
To quantify these statements, let us for a moment con sider the 0{p^) HBChPT calculation. If one switches off the 



complete contribution from the axial correlator, Eq.(4.4), the singlet rate is enhanced by a factor of 1.7, whereas the 
triplet is decreased by a factor of about 170 ! Consequently, we then have Ftot — 7 x 10"^ s"^, which is an order of 
magnitude smaller than the value given in the previous section. If one on the other hand sets = k„ = in the V, 
VV and AV correlators, the singlet rate increases by a factor of about 2.7 and the triplet rate decreases by a factor 
of 1.3, leading to a total rate of 68xl0-3s-^. 

It is also instructive to consider the chiral limit, m^r = 0. To be specific, we discuss the C'(p^) HBChPT calculation, 
with the C(e^) SSE results being very similar. To leading order in l/M^v, one encounters a pole at a; = m^/2 in the 
A and AV correlators. This can be seen by looking at pion pole terms in the chiral limit, which take the form 

^ ^ (4.46) 



(/ - ml - 2mf,l'o 



and using the leading order result Iq — — u + (^(l/Mjv), cf. Eq.(4.38). This well-known Bethe-Heitler pole is 
independent of the angular variable x. The condition for this pole to appear is that the pion mass has to be below 
the muon mass. For zero pion mass, there is another pole at the same energy for x — —1 stemming from the terms 
~ (r — r')~^. Furthermore, the chiral limit photon spectra of the 0{p^) HBChPT calculation are shown in Fig.^. We 
remark that these spectra are very different from the ones with the physical pion mass due to the abovementioned 
singularities. We also point out that the pion mass effects are larger in RMC than in the OMC case discussed above. 
In particular, setting e.g. ttIt^ = 125 MeV, the rate in the para state (total rate) increases to 99xl0-'^s-^. We will 
discuss the implications on the TRIUMF measurement in the following section. 
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3. Discussion of the TRIUMF result for gp 



The photon spectra discussed in section IV D 2 aUow in principle to determine the induced pseudoscalar form factor. 
The TRIUMF result for gp '\s obtained by multiplying the terms proportional to the pseudoscalar form factor with a 
constant denoted R (the momentum dependence assumed to be entirely given by the pion pole). The value of R is 
then extracted using the model of Fearing et al. to match the partial rate for photon energies larger than 60 MeV. 
If we perform such a procedure, we get a similar shift in the partial photon spectra (using the same weight factors 
for the various n — p states as given in ref. Q). It is, however, obvious from our analysis that such a procedure is 
not legitimate. By artific ially enhancing the contribution ~ gp (to simulate this procedure, we have introduced the 



factor R in Eqs.(4.4,4.12)), one mocks up a whole class of new contact and other terms not present in the Born term 
model. To demonstrate these points in a more quantitative fashion, we show in fig.^ the partial branching fraction 
for our calculation in comparison to the one with gp enhanced by a factor 1.5 and a third curve, which is obtained 
by increasing gp only by 15%, take hi — 24 (i.e. enhancing this couphng by a factor of two) and use A = 273 MeV, 
since in the dispersion theoretical analysis of pion-nucleon scattering the pole in the P33 partial wave is located at 
W = 1210 MeV. This is shown in fig.§ by the dashed line and it shows that such a combination of small effects can 
explain most (but not all) of the shift in the spectrum. This is further sharpened by using now the neutral pion mass 
of 134.97 MeV instead of the charged pion mass, leading to the dotted curve in fig.^. Since the pion mass difference 
is almost entirely of electromagnetic origin, one might speculate that isospin-breaking effects should not be neglected 
(as done here and all other existing calculations). Furthermore, as discussed above a slight change in the occupation 
numbers pi/2 and P3/2 would also lead to an increase in the ortho capture rate which could close the gap between 
the empirical and theoretical results. For example the dashed curve in fig.|| would be moved from 0.42 to 0.48 while 
the dotted one would go from 0.50 to 0.54 with pi/2 — 0.95 and — 0.05. The situation is reminiscent of the 
sigma term analysis, where many small effects combine to give the sizeable difference between the sigma term at zero 
momentum transfer and at the Cheng-Dashen point. We further point out that although the present investigation 
combined with the findings in ref. [|6| does not seem to give any large new term at order or from the chiral loops 
at third order, it can not be excluded that one-loop graphs with insertions from the dimension two chiral Lagrangian 
(which are formally of fourth order) can give rise to larger corrections than the third order loop and tree terms 
calculated in ref. ||l6[. In fact, as we noted before, the photon spectrum is more sensitive to changes in the anomalous 
magnetic moments than to the induced pseudoscalar coupling. Therefore it can be speculated that one-loop graphs 
with exactly one insertion ~ k„ can generate large corrections. This appears plausible but needs to be supported by 
a fourth order calculation. Such a mechanism would, however, be much more natural than the simple rescaling of 
gp based on tree level diagrams only. Another point against this simple rescaling comes from OMC. Indeed if this 
rescaling holds for RMC it should also hold for OMC. We thus have performed a similar calculation in OMC. Taking 
the same value for R, one would obtain 

rg^^c' - 172.8 X .-1, r^^-M^c''^^ ' = 634.6 X ,s-\ T'^^^='-' = x s-\ (4.47) 

leading to r°^°' = 477 x s^^, which is lower than the error bars of the experimental result from Bardin et al. 

. As expected, the singlet and triplet capture rates are much more sensitive to the details of the interaction than 
the total rate. As a conclusion we note that the effect of enhancing the capture rates in RMC via setting i? = 1.5 
leads to a strong reduction of the corresponding OMC rates leading to conflicts with the experimentally determined 
ortho capture rate. 

To summarize this discsussion, we have pointed out that two effects in particular have to be investigated in more 
detail: 

- the occupation numbers of the atomic structure in muonic atoms/molecules need to be carefully re-examined 
by experts in this field. 

- the N^LO calculation should be redone including all isospin breaking effects because of the sensitivity to the 
exact pion mass in the pion-pole contributions, for example. 

The sum of these small effects should explain the observed photon spectrum, as we believe that the proper 
hadronic/weak physics part is well under control by now, as our analysis has re-confirmed. A simple rescaling of 
the pseudoscalar coupling constant should no longer be considered. 

V. SUMMARY 

In this manuscript, we have considered ordinary and radiative muon capture on the proton in the framework of 
the small scale expansion to third and second order in small momenta, O(e^) and O(e^), respectively. We have also 
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discussed the induced pseudoscalar form factor of the nucleon and its determination from the TRIUMF RMC data. 
The pertinent results of this investigation can be summarized as follows: 



(i) To third order in the small scale expansion, ordinary muon capture is almost not affected by A(1232) isobar 
effects. The only effect comes via the Pauli radius and is extremely small. The NLO contribution to the total 
capture rate amounts to a 25% correction of the leading term. This is in agreement with naive dimensional 
counting, which lets one expect corrections of the size m^/K^. This calculation involves very few and well- 



controlled parameters. We argued that the formula derived from a relativistic calculation, see Eq.(4.27), and 
extensively used in the literature does not hold in the modern view point of power counting. We have stressed 
the importance of the upcoming PSI experiment p^ . 

(ii) To second order in the small scale expansion, we have considered radiative muon capture. A(1232) related 
effects only appear at NLO and its effects on the total capture rate and the photon spectrum are of the 
order of a few per cent. T he smallness of the A(1232) contributions is due to a combination of effects as 
discussed in section IV D 2 . This agrees with earlier findings in a more phenomenological approach . Isobar 



effects can therefore not resolve the discrepancy between the TRIUMF measurement for the partial decay width 
r(u; > 60MeV) and the theoretical predictions. We have, however, pointed out severe loopholes concerning 
the extraction of gp as done in ref. one of them being the contradiction with the OMC data. In our 
opinion the most probable explanation of the discrepancy is a combination of many small effects, as detailed in 



section IV D 3 



(iii) The induced pseudoscalar form factor measured in charged pion electroproduction is not very well determined, 
but clearly is in agreement with the one-loop chiral perturbation theory prediction Q. A more precise mea- 
surement for small invariant momentum transfer squared is called for |3q]. 
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APPENDIX A: FEYNMAN RULES 

In this appendix, we collect the pertinent Feynman rules for calculating OMC and RMC. These read: 
Isovector vector source in (g^), nucleon: 

« ■ V + ^ (in + rf) ■^~in + rf)■vv■^r) + ^J^ + k,)[S■^^,S■q] , (A.l) 



Isovector axial source in, nucleon: 



i2gAS-a-i^S-{n+rf)va, (A.2) 



Isovector axial source in, pion out (fc^): 



Fo i Tr[r' • fc + a • fc r*] , (A.3) 



Isovector axial source, vector source, pion: 



2 vector sources in, nucleon: 



2Mn 



^Tr{[v^,r^]a'^ +a^[v^ri} , (A.4) 



{v + v('^)'~\viv + v^'y)Y -il + k,)[S^,S''][v,,v,] \ , (A.5) 
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vector source and isovector axial source in, nucleon: 



gA 
Mo 



5 • (v + v'^'^A V & + V aS ■ 



in, nucleon and vector source (fc^) out: 



bi 



(A.6) 



(A.7) 



Nucleon in. A' and vector source (fc^) out: 



. 61 



2Mf 



Tr[T*(fc'V'^-rv'^)]S'^ , 



(A., 



Isovector axial source, nucleon. A' : 



«57rArATr[T*ap] . 



(A.9) 
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Tables 



Quantity 


Symbol 


Value 


Units 


Proton mass 




938.27321 


MeV 


Neutron mass 


M„ 


939.56563 


MeV 


Nucleon mass 


Mn 


938.91942 


MeV 


Pion mass 




139.56995 


MeV 


Muon mass 




105.658389 


MeV 


Axial-vector coupling 




1.2670 




Pion-Nucleon coupling constant 


g-irNN 


13.05 




Pion decay constant 




92.42 


MeV 


isovector Dirac radius 


rl 


(0.585)1/2 


fm 


0{t') SSE isovector Pauli radius 


r5(SSE) 


(0.61)1/2 


fm 


C)(p^)HBChPT isovector Pauli radius 


r5(HBChPT) 


(0.52)1/2 


fm 


isovector axial radius 


VA 


0.65 


fm 


Fine-structure constant 


a 


1/137.0359895 




Fermi constant 


Gf 


1.16639 • 10"^ 


GeV-2 


CKM matrix element 


Vud 


0.9740 





TABLE I. Values of the various masses, couplings and other constants used in the text. 
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Figures 



V 



'xz ^ 



FIG. 1. Graphical representation of the V, A (first row), and VV (second row) correlators. Nucleons, pions, photons and 
external gauge fields are denoted by solid, dashed, wiggly and zigzag lines, in order. 
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FIG. 2. Nucleon contributions to the VA correlator. All external gauge fields are axials. 
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FIG. 3. Delta contributions to the VA correlator. All external gauge fields are axials. 




FIG. 4. 0{e^) SSE photon spectra for RMC for the singlet, triplet, para (statistical) and ortho states of the fJ. — p system. 
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FIG. 5. Relative change in the photon spectra for RMC due to expUcit A(1232) related effects for the singlet, triplet, para 
(statistical) and ortho states of the n — p system. 




FIG. 6. 0{p^) HBChPT photon spectra for RMC for the singlet, triplet, para (statistical) and ortho states of the n — p 
system in the chiral limit. 
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FIG. 7. The "world data" for the induced pseudoscalar form factor Gp{t). Dashed curve: Pion-pole prediction. Solid curve: 
0{e^)/0{jp^) SSE/HBChPT prediction. The pion electroproduction data (crosses) are from ref. M. Also shown is the OMC 
result a± t — — 0.88m^j from ref. (open square). 




FIG. 8. Photon spectra for RMC for the branching ratios of the singlet, ortho and para states as used in the TRIUMF 
analysis. Solid line: Prediction of the small scale expansion to order e^. Dash-dotted line: Same as the solid line but with gp 
scaled by a factor 1.5. Dashed line: Various small modifications as explained in the text. Dotted line: Same as the dashed line 
but using the neutral instead of the charged pion mass. 
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